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ABSTRACT 

We present a detailed analytical study of the propagation of ultra high energy (UHE) particles in 
extragalactic magnetic fields. The crucial parameter which affects the diffuse spectrum is the separation 
between sources. In the case of a uniform distribution of sources with a separation between them 
. much smaller than all characteristic propagation lengths, the diffuse spectrum of UHE particles has a 

universal form, independent of the mode of propagation. This statement has a status of theorem. The 
proof is obtained using the particle number conservation during propagation, and also using the kinetic 
equation for the propagation of UHE particles. This theorem can be also proved with the help of the 
Oh' diffusion equation. In particular, it is shown numerically, how the diffuse fluxes converge to this universal 

, spectrum, when the separation between sources diminishes. We study also the analytic solution of the 

diffusion equation in weak and strong magnetic fields with energy losses taken into account. In the 
Q^ \ case of strong magnetic fields and for a separation between sources large enough, the GZK cutoff can 

practically disappear, as it has been found early in numerical simulations. In practice, however, the 
source luminosities required are too large for this possibility. 

ly-^ ' Subject headings: UHE Cosmic rays, propagation of cosmic rays, GZK cutoff. 

a^ : 

O ' 1. INTRODUCTION 

■ 

■ The propagation of UHECR protons and nuclei with i? > 1 x 10^^ eV in the large-scale intergalactic magnetic field 
, (IMF) remains an open problem, mainly because the knowledge of the IMF is still very poor. The possibilities vary 

■ between rectilinear propagation in a weak field and diffusive propagation in a strong magnetic field. The experimental 
\ data on IMF and the models of origin of these fields do not allow at present to choose even between the two extreme 

O ^1 propagation regimes mentioned above. 

I ■ Most reliable observations of the intergalactic magnetic field are based on the Faraday rotation of the polarized radio 
O \ emission (for the reviews see Kromberg (1994), Valle (1997), Carilli and Taylor (2002)) . The upper limit on the Faraday 
> • rotation measure (RM) in the extragalactic magnetic field, obtained from the observations of distant quasars, gives an 
, upper limit of RM < 5 rad/m . It implies an upper limit on the extragalactic magnetic field on each assumed scale of 
• • • coherence length (Kromberg (1994), Valle (1997), Ryu et al. (1998)). For example, according to Blasi et al. (1999a) for 
an inhomogeneous universe Bi^ < 4 nC on a scale of coherence Ic = 50 Mpc. 

According to observations of the Faraday rotations the extragalactic magnetic field is strongest, or order of 1 fiG, in 
^ clusters of galaxies and radiolobes of radiogalaxies (Valle (1997), Kromberg (1994), Carilli and Taylor (2002)). The largest 
Ci scale in both structures reaches Zc ~ 1 Mpc. Most probably various structures of the universe differ dramatically by 
magnetic fields, with very weak field in voids and much stronger in the filaments (Ryu et al. (1998)). Superclusters seem 
to be too young for the regular magnetic field to be formed in these structures on a large scale Ic ^ 10 Mpc. 

In case of hierarchical magnetic field structures in the universe, UHE protons with E > Ax 10^^ eV can propagate in a 
quasi-rectilinear regime. Scattering of UHE protons occurs mostly in galaxy clusters, radiolobes and filaments. Defiections 
of UHE protons can be large for some directions and small for the others. The universe looks like a leaky, worm-holed 
box, and correlation with the sources can be observable (see Tinyakov and Tkachev (2001)), where correlations of UHECR 
with BLLacs are found). Such a picture has been suggested by Berezinsky et al. (2002b). 

A promising theoretical tool to predict the IMF in large scale structures is given by magneto-hydrodynamic (MHD) 
simulations. The main uncertainty in these simulations is related to the assumptions concerning the seed magnetic field. 

The MHD simulations of Sigl et al. (2004) and Sigl et al. (2003) favor the hierarchical structure with strong magnetic 
fields. Assuming an inhomogeneous seed magnetic field generated by cosmic shocks through the Biermann battery 
mechanism, the authors obtain a ~ 100 nG magnetic field in filaments and ~ 1 nG in voids. In some cases they consider 
IMF up to a few micro Gauss as allowed. In these simulations UHECR are characterized by large deflection angles, of 
the order of 20°, at energies up to i? 10^*^ eV (Sigl et al. (2003), Sigl et al. (2004)). Thus, the scenario that emerges in 
these simulations seems to exclude the UHECR astronomy. These simulations have some ambiguity related to the choice 
of magnetic field at the position of the observer (Sigl et al. (2003), Sigl et al. (2004)). The authors consider two cases: a 
strong local magnetic field B ~ 100 nG and a weak field B ^ 100 nG. The different assumptions about the local magnetic 
field strongly affects the conclusions about UHECR spectrum and anisotropy. 

The essential step forward in MHD simulations has been made recently. In Dolag et al. (2003) the Local Universe is 
simulated with the observed density and velocity field. This eliminates the ambiguity for the local magnetic field, that 
is found to be weak. The seed magnetic field, used in this simulation, is normalized by the observed magnetic field in 
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rich clusters of galaxies. The results of these constrained simulations indicate the weak magnetic fields in the universe 
of the order of 0.1 nG in typical filaments and of 0.01 nG in voids. The strong large-scale magnetic field, B ~ 10'^ nG, 
exists in clusters of galaxies, which, however, occupy insignificant volume of the universe. The picture that emerges from 
simulations of Dolag at el. (2003) favors a hierarchical magnetic field structure characterized by weak magnetic fields. 
UHE protons with E > Ax 10^^ cV can propagate in a quasi-rectilinear regime, with the expected deflection angles being 
very small < 1°. However, until direct observational evidences for this picture becomes available, an alternative case of 
propagation in strong magnetic fields, with diffusion as an extreme possibility, can be hardly excluded. 

This case has been studied in Sigl et al. (1999), Lemoine et al. (1999), Stanev (2000), Harari et al. (2002), Yoshiguchi 
et al. (2003), Deligny et al. (2003). An interesting features found in these calculations are small-angle clustering of UHE 
particles as observed by Hayashida ct al. (1996), Hayasliida et al. (1999), Uchiori ct al. (2000), Glushkov and Pravdin 
(2001), and absence of the GZK cutoff in the diffusive propagation, when the magnetic field is very strong. Many aspects 
of diffusion of UHECR have been studied in numerical simulation by Casse et al. (2002). 

We shall illustrate UHECR propagation in strong magnetic fields by the calculations by Yoshiguchi et al. (2003). 
The authors performed MC simulations for propagation in random magnetic field with the Kolmogorov spectrum of 
turbulent energy density Wk oc fc", with n = —5/3. The basic (largest) coherent scale is chosen as Ic with the mean field 
Bq. Numerically these parameters vary in the range (1 - 40) Mpc for Ic, and (1 - 100) nG - for Bq. The sources are 
taken as galaxies from Optical Redshift Survey catalog with absolute magnitude Mb brighter than some critical value 
Mc- The calculated quantities are the energy spectrum, anisotropy and small-angle clustering. The observed small-angle 
clustering and absence of the GZK cutoff in the AGASA observations can be reproduced in the case of strong magnetic 
field B > 10 nG. 

Diffusive propagation of extragalactic UHECR has been studied already in earlier work. The stationary diffusion from 
Virgo cluster was considered by Wdowczyk and Wolfendale (1979), Giller et al. (1980) and non-stationary diffusion from 
a nearby source was studied by Bcrczinsky ct al. (1990a), Blasi and Olinto (1999b) using the the Syrovatsky solution 
(Syrovatskii (1959)) of the diffusion equation. In this case the GZK cutoff can be absent. A very similar problem was 
considered again more recently by Isola et al. (2002)). 

In this paper we shall study how propagation influences the diffuse energy spectrum of UHECR. We shall prove the 
theorem that if distance between sources is much smaller than all propagation lengths, the spectrum has the same 
universal form independent of the mode of propagation. For diffiision in magnetic flelds we shall demonstrate how the 
spectra converge to the universal one when the separation between sources diminishes. Finally, we shall obtain, with the 
help of the Syrovatsky solution, the spectra for strong magnetic field. In this case with large enough separation between 
sources, the GZK cutoff becomes weak or absent, and we will discuss the physical explanation of this phenomenon. 

2. PROPAGATION THEOREM 

Let us consider a case when identical UHECR sources are distributed uniformly ^ in the space, with d being the 
separation between sources. We will demonstrate that if d is less than all other characteristic lengths of propagation, such 
as diffusion length ld{E) and energy attenuation length latt{E) given by 

latt = cE/idE/dt), (1) 

then the diffuse energy spectrum has a universal (standard) form independent of the mode of particle propagation. In 
particular, under the conditions specified above the magnetic field, both weak and strong, does not affect the shape of 

the energy spectriim. 

Explicitly, this theorem can be formulated as follows: 

For a uniform distribution of identical sources with separation much less than the characteristic propagation lengths, 
the diffuse spectrum of UHECR has a universal (standard) form, independent of the mode of propagation. 

First, we shall consider the proof based on the conservation of the number of particles (e.g. protons or nuclei) during the 
propagation. Let t be the age of the universe with the present age taken as to- The number of particles per unit volume 
of the present universe is equal to the number of particles injected into this volume during all history of the universe, 
independent of the mode of propagation. The homogeneity of particles needed for this statement is provided by almost 
homogeneous distribution of sources. Thus, the comoving space density of particles np{E) from uniformly distributed 
sources with an age-dependent comoving density ns{t) and age-dependent generation rate by a source, Q{E,t), is given 

by 

np{E)dE= [ " dtQ{Eg,t)ns{t)dEg, (2) 

where Eg{E,t) is the required generation energy at age t, if the observed energy is E. Eq. (2) does not depend on the 
way particles propagate. 

The homogeneous distribution of particles in presence of inhomogeneous magnetic fields follows from the Liouville 
theorem and can be explained in the following way. Suppose we have an observer in the space with a magnetic field. 
The diffuse flux in any direction is given by the integral / ns{l)dl over the trajectory of a particle (or antiparticle emitted 
from the observation point). If is almost homogeneous, the integral depends on the time of propagation and does not 
depend on the strength and inhomogeneity of magnetic field. 

^ In this paper we shall distinguish between uniform and homogeneous distribution of the sources: under the latter we assume a continuous 
and distance-independent distribution. 
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To find the explicit form of the universal spectrum np{E), one needs some additional assumptions. Let us consider the 
protons as primaries with continuous energy losses due to interaction with the CMB radiation 

dE/dt = -b{E,t). (3) 

Here and everywhere below b{E, t) > 0. We shall use the connection between the redshift z and cosmological time t 
according to the standard cosmology 

dt = , (4) 

Ho{i + z)^{i + zra^ + nA 

with Hq, rj„, and Qa being the Hubble constant, relative cosmological density of matter and relative density of vacuum 
energy, respectively. 

The generation rate is assumed to be the same for all sources and is taken in the form 

Q{E, t) = Lp{l + zrKijg)qg,,,{Eg) (5) 

where Lp is the CR luminosity of the source with (1 + z)" describing the possible cosmological evolution of luminosity. 
The normalization factor K{'yg) is 7g — 2 if 7g > 2 and 1/ ln(iJmax/-Bo) if 7g = 2, with Eq and E'max the minimum and 

maximum generation energies, respectively. The comoving density ns{t) of the sources can also contain the evolutionary 
factor (1 + z)'^, with a + f3 — m. Here and everywhere else we assume Eq = 1 GeV, all energies E are measured in GeV 
and Lp in GcV/s. 

Then from Eq. (2) we obtain the normalized universal spectrum as: 

Jp{E) = ^CoK{^g) dz (1 + zTqgUE,)'^ (6) 

where Jp is diffuse flux, £o = LpUs is the emissivity at z = 0, dt/dz is given by Eq. (4), Eg is calculated as Eg = Eg{E, z) 
and dEg/dE is given by (Berezinsky and Grigorieva (1998), Berezinsky et al. (2002a)) 

'^•^'•>-i.^..)e.,{^r,. .„ (g#i) ]. m 



where h^iE) is the proton energy loss at z = 0. 

The generation spectrum qgen{Eg) is not known apriori. Three kinds of spectrum might be considered: (i) the power 

law spectrum qgen{Eg) = Eg with 7g > 2 and the normalization factor K{jg) =7^ — 2, (ii) the power-law spectrum 
qgen{Eg) = 1/Eg with K(pig) = 1 / \n{E„iax / Eo) and (iii) the complex spectrum 

fenl^gj - I E-\Eg/E,)-^9 at Eg > E, 

with 

^(7.) = T-ET^ 



In numerical calculations with the complex spectrum (8) we use 7g = 2.7, which gives the best fit to observational data 

(Berezinsky ct al. (2003)). 

The universal spectrum (6) with ggon from Eq. (8) and with m = is shown in Fig. 1 in comparison with experimental 
data of AGASA and HiRes arrays. From these figures one can see the good agreement of the universal spectrum with 
experimental data of both detectors up to 8 x 10^^ eV. The required emissivities jCq are: £o — 1.8 x lO"'^ erg/Mpc^yr and 
Co 8.9 X 10^^ erg/Mpc^yr for the AGASA and HiRes data, respectively. The excess of AGASA events at £ > 1 x lO^o eV 
needs another component of UHECR. 

Under the condition d << /prop the energy spectrum (6) is the same (universal) for rectilinear propagation and propa- 
gation in weak and strong magnetic fields. If /prop < d (or /prop « d) the propagation theorem is not valid any more. The 
distortion of the universal spectrum occurs at E > Eqzk due to several effects: inhomogeneity of distribution of UHECR 
sources (local enhancement or local absence of the sources), fluctuation in the distribution of sources, fluctuations of the 
photopion energy losses over the distance d and due to the small diffusion length < in case of strong magnetic fields. 

3. KINETIC EQUATION AND UNIVERSAL SPECTRUM 

We have obtained the universal spectrum above from the very general condition of particles conservation during prop- 
agation. In this Section we shall demonstrate that the universal spectrum follows also from the kinetic equation. 

Let us consider the kinetic equation for a homogeneous distribution of the sources (see footnote ^). Since in this case 
the particle distribution is also homogeneous, the diffusion term in the kinetic equation is absent, and it has the form: 

^^^4r^ = ^ [b{E, t)np{E, t)] + Qg{E, t), (9) 
with Qg = UsQ. Let us start with a simple stationary equation when Qg, Up and b do not depend on time: 

-^[h{E)np{E)]=Qg{E) (10) 
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The solution of Eq. (10) is 



(11) 



where Eg{E, t) is the generation energy at time t and we used dEg = —h{Eg)dt. Using dE = —b{E)dt for the same interval 
dt we obtain b{Eg)/b{E) = dEg/dE and 

r ^TP 

(12) 



n,{E)= I dtQg{Eg{t))^, 



in agreement with Eqs. (2) and (6). 

Let us now come back to Eq. (9) with time dependent quantities Qg,np and &, where b{t) includes also adiabatic 
energy losses due to redshift. Introducing the new quantities: b{E,t) = b{E)/H{t), Qg{E,t) = Qg{E,t)/H{t) and 
T = (1 + Zmax)/{^ + z), where Zmax is the maximal redshift in the evolution of sources, we obtain the equation 



dnp{E,T) _ db{E,T) 



dr dE 
Eq. (13) may be solved with the help of an integration factor 



np{E,T)+Qj,{E,T). 



,^r)=e.J-[dr'§ir') 



Assuming np{E, Zm) = 0, we obtain 



np{E,T) 



m(t) 



r dT'n{T')Qg{E,T') 

.Jo 



(13) 



(14) 



(15) 



At this stage the energy losses should be separated into those due to redshift and those due to the interaction with the 
CMB radiation 

b{E', z) = E'HoVnmil + z)^ + Oa + (1 + zfb^^"" {E'{1 + z)) (16) 

where E' is an arbitrary energy at epoch z. 

Coming back to the variable z in Eq. (13) and inserting there db{E, z)/dE from Eq. (16) we obtain after integration 



X exp 



- r 

-f^o Jo 



Vf2„(l + 2)3 + QA 

dz'{l + z'f (dbo{E',z') 



^^m{l+z'f + ^K \ dE' 

/ E'={l+z)Eg{z') 



(17) 



which, after introducing the explicit expression for Qg from Eq. (5), coincides exactly with Eqs. (6) and (7). 
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4. DIFFUSION EQUATION AND THE SYROVASTSKY SOLUTION 

In the previous section we have considered the case of a homogeneous distribution of sources, and proved that the 
universal spectrum is a solution of the kinetic equation (9). Wc shall now assume a non-homogeneous distribution of 
the sources. Then one should add the propagation term in Eq. (9). We assume the diffusive propagation in large-scale 
magnetic fields. For a single source at point fg the diffusion equation has the form: 

dnp{E,r) ,. ,^ „ d 



■div [D{E,r, t)\Inp{E,\ 



[b{E)np{E,r)]=Q{E,t)5{r^rg 



(18) 



where D{E, r, t) is the diffusion coefficient that depends on the magnetic field structure. We refer the reader to the Section 
5 for a detailed discussion of the relation between the diffusion coefficient and the magnetic field. 

In the case when D, b and Q depend only on energy, the exact analytic solution of Eq. (18), found by Syrovatskii (1959), 

is 

> exp 



ip{E,r) 



1 



/ 



"4A(E,Eg) 



where 



X{E,E,) 



f 

Je 



de 



{inX{E,Eg)) 
6(e) 



3/2 ■ 



(19) 



(20) 



is the Syrovatsky variable which has the meaning of the squared distance traversed by a proton in the observer direction, 
while its energy diminishes from Eg to E. 

According to the propagation theorem, integrating the Syrovatsky solution over homogeneously distributed sources with 
density rig. 



m 



f 



dEgQiEg) 



exp 



'iHE,E,) 



{A'KX{E,Eg)) 



3/2 



(21) 



must result in the universal spectrum. One easily sees this by changing the order of integration in Eq. (21) and using 



f 

Jo 



dr- 



(47rA) 



3/2 



exp 



r 

4A 



which gives the diffuse space density of protons 

np{E) 



Us r 



dEgQ{Eg) 



(22) 



b{E) JE 

in agreement with Eq. (11), where Qg{E) = nsQ{E) 

Coming back to the non- homogeneous case, let us consider a lattice with total size a and UHECR sources located in 

the lattice vertexes separated by a distance d. Using this model for the source distribution, we shall demonstrate the 
consistency of the Syrovatsky solution with the propagation theorem. We will show that when the separation between 
sources, d, tends to a small value dum, the Syrovatsky solution gives the universal spectrum. 
Using the lattice distributed sources and the Syrovatsky solution, one obtains the diffuse flux as 



Jp{E) = 



1 



47rb{E)^JE 



E 



dEgQiEg)- 







exp 


4\{E.Eg)_ 



(23) 



{47rX{E,Eg)f' 

where X{E, Eg) is given by Eq. (20) and the summation goes over all sources in the lattice vertexes. The maximum energy 
in Eq.(23) is given by the smaller of two quantities: the maximum acceleration energy E™^^, and the generation energy 
Eg{E,cto) of a proton with present energy E propagating during a time tQ. For the energies E at interest, the former is 
always smaller that the latter, and we use as i^/max the maximum acceleration energy. 

We have to specify now the diffusion coefficient D{E), which determines X{E, Eg) in Eq. (23). Putting off the detailed 
discussion until the Section 5, we will give here a short description of the diffusion coefficients used in this work. 

We assume diffusion in a random magnetic field with the mean value Bq on the maximum coherent length Ic- This 
assumption determines the diffusion coefficient D{E) at the highest energies when the proton Larmor radius, rL{E) » Ic- 



■HE) 



D{E) = 

At "low" energies, when rL{E) < 1^ we shall consider three cases: 

(i) energy-independent diffusion coefficient 

D = gC/c, 

(ii) the Bohm diffusion coefficient, which provides the lowest value of D 

Db{E) = \crL{E), 



(24) 



(25) 



(26) 
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Fig. 2. — Convergence of the diffusion spectrum to the universal spectrum in the case of the Bohm diffusion, D = Dq{E/Ec) at E < Ec 
(right panel) and in the case of diffusion with D = Dq at B < Be (left panel). In the calculations the complex generation spectrum (8) has 
been used. The separations d of the sources are indicated in the figure. The spectra are compared with the AGASA data. 



(iii) the Kolmogorov diffusion coefScient 

In all three cases we normalized the diffusion coefficients by (l/3)c^c at tl = h (see Section 5). The characteristic 
energy Ec of the transition between the high energy and low energy regimes is determined by the condition r^iE) = 1^ 
and is 

One can describe the low-energy and high-energy diffusion regimes with the help of an interpolation formula for the 
diffusion length: 

14E) =Ad + (29) 

with Arf = Ic for the regime with D = const, = rL{E) for the Bohm diffusion and A,^ = Ic{tl/IcY^^ for the Kolmogorov 
regime. For the later use, we shall formalize the description of these three regimes using 

Ad(S) = Utl/IcT, (30) 
with a equal to 0, 1 and 1/3 for the D = const, Bohm and Kolmogorov regimes. 
For completeness we shall give also the numerical expression for the Larmor radius: 

r^E) = 1.08 X 102__|^1|^ Mpc. (31) 

One can see, therefore, the existence of two different propagation lengths that should be compared with the distance 
d between sources. These two lengths are latt given by Eq. (1) and Id given by Eq. (29). The former is large even at 

the highest energies (latt ^ 25 Mpc at E ^ 10'^^ eV) while the latter can be small enough at the energies of interest 
10^^ < E < 10^^ eV. For a representative case Bq — 100 nG and ~ I Mpc we have ld{Ec) = ?c = 1 Mpc at 
Ec ~ 1 X 1020 eV. 

In Fig. 2 we compare the universal spectrum with diffusion spectra characterized by different separations d. The 
diffusion spectra are computed using the diffusion length (29) for the two extreme cases A^ = rL{E) and A^ = Ic, both 
for the representative ease Bq = 100 nG and Ic = 1 Mpc. Here and below we use for the total lattice size a = 300 Mpc, 
while the control calculations have been performed up to a = 1000 Mpc. Increasing a does not change the fluxes. The 
universal spectrum is calculated with time-independent energy losses, as in the Syrovatsky solution, according to Eq. (6) 
with Co = Lpns 2.8 x lO-^*^ erg/Mpc^yr. 

One can see that in both cases the diffusion spectrum tends to the universal one as the separation d between sources 
diminishes. 

5. DIFFUSIVE PROPAGATION 

In this section we will study the propagation of UHE particles in the intergalactic turbulent magnetic plasma. Turbulent 
motion can be considered as random pulsations at different scales, described as an ensemble of many waves with large 
amplitudes and random phases. The energy is assumed to be injected at the largest scale Ic due to the external processes, 
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such as the contraction of large-scale structures, large-scale shocks etc, and the scale Ic might be determined by these 
processes. This scale might coincide with the size of the region filled by the turbulent plasma, in case of small structures. 
The smallest scale Zmin is determined by the dissipation of turbulent motion into thermal energy. 

We shall consider magnetic turbulence described by the superposition of MHD waves >lexp[i(— -|- kr)] with different 
frequencies u) and different wave numbers k. The amplitude A is related to magnetic field B, density p, pressure P, 
velocity v etc. These waves propagate with the Alfven velocity ua = B/y/Airp. The kinetic energy of gas in turbulent 
motion and the magnetic energy in MHD waves are equal at all scales A (Landau and Lifshitz (1987)): 

pvl ~ Bl/^TT. (32) 

For a detailed description of magnetic turbulent plasma and its description as an ensemble of MHD waves we refer the 
reader to the books: Landau and Lifshitz (1987), Arzimovich and Sagdeev (1979). 

For practical calculations of the UHECR propagation, the magnetic field scale spectrum {B^)k as a function of k is 
needed. One should clearly distinguish between the magnetic field B\ = {B)x, on the scale A = 2n/k and the Fourier 
component of the field B{k) (see Appendix). The physical magnetic field which determines the diffusion of UHE particles 
is Bx. 

The turbulence spectrum is described by the dependence of the spectral energy density w{k) on the wave number k. 
It can be expressed using the magnetic energy density kw{k) ~ Bf/Sw. Usually the spectra have a power-law form: 

w{k) oc 

In many practical cases the observations confirm the Kolmogorov (Kolmogorov (1941)) spectrum of turbulence with 
m = 5/3. This spectrum belongs to a class where the initial energy is assumed to be injected at the largest scale Ic and 
then is transferred to lower scales I in non-linear processes of wave interactions until it dissipates to thermal energy at 
the lowest scale /min- The Kolmogorov spectrum can be derived assuming that energy flux from one scale to another does 
not depend on the scale k (Arzimovich and Sagdeev (1979)), i.e. 

kw{k)/Tk = const, (33) 

where is the time for energy transfer to the scale k. Since this process is caused by the non-linear term {vV)v kv"^ 
in the Euler equation, one can estimate Tk as 

Tfe ~ 1/kvk, (34) 

where the turbulent velocity can be estimated from equation kw{k) ~ pv^. Then one immediately obtains w{k) oc k~^^^, 
i.e. the Kolmogorov spectrum. 

Landau and Lifshitz (1987) argue that while scale-independent energy flux is a reasonable assumption in ordinary 
hydrodynamics, for MHD waves it is proportional to vf.. In this case one obtains the Kraichnan (1965) spectrum w{k) oc 

In the case of shock waves, the spectrum is (Vainstein et al (1989)) w{k) oc k^. 

5.1. Diffusion coefficient for UHE protons 

Now we shall proceed with calculations of the diffusion coefficient D(E) for UHE protons in extragalactic magnetic 
fields, which are characterized by the spectrum of turbulence w{k) oc fc""*, with the basic scale Ic and with the magnetic 
field on this scale Bq. We shall use the characteristic energy Ec, defined by Eq. (28) from the condition rL{Ec) = Ic, 
where the Larmor radius is given by Eq. (31). 

First we obtain the asymptotic expression for D{E) valid when r^ ^ Ic, i.e. at E ^ Ec- The diffusion length is 
determined as the distance at which the average angle of scattering satisfies 

r)~„(i^)^i. 

where the number of scatterings n is given by n ^ Id/lc- It results in the diffusion length ld{E) ~ r^/lc, and thus the 
diffusion coefficient, D{E) = {l/3)clci{E), in the asymptotic limit tl ^ Ic is given by 

D_,(£;) = 3.6 X 10-i?|o (laZ^) ' (i^) cmVs. (35) 

We shall calculate now the diffusion coefficient in the low-energy limit <C Ic (i-e- a.t E <^ Ec). We will follow the 
V. Ptuskin method given in Chapter 9 of the book Berezinsky et al. (1990b). Let us consider the scattering of UHE 

protons with ^ in a magnetic field of MHD wave B\exp[i{kf — u)t)]. The magnetic field on the basic scale Bq is 
considered as constant field for smaller scales A. 

Scattering of UHE protons off the MHD waves in the regime r^ <SC Ic is dominated by resonance scattering (see Lifshitz 
and Pitaevskii (2001), Sections 55 and 61). The condition for resonant scattering is given by u)' = scu'q, where iv' and 
u}'^ = eB/^'mc are the wave frequency and the giro-frequency in the system K' at rest with particle motion along the 
field Bq, and s = 0, ±1, ± 2.... After a Lorentz transformation to the laboratory system, one obtains 

oj — k\\Vz = sivb, (36) 
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where lu is the wave frequency, fcy is the projection of the wave vector onto the direction of Bq , Vz is the projection 
of particle velocity onto the same direction, and ojb = eB/^mc. For a magnetized plasma ri, <i; A_l (or equivalently 
k±v±/uiB 1) and s = ±1 (Berezinsky et al. (1990b)). Thus from Eq. (36) one derives the resonant wave number 



Lo{k) ± OJb 



^ = —, (37) 



Vz 

where fj, = cos6 and vl = v/cob is the Larmor radius. In deriving Eq. (37) we have used u){k) <C cjb, which follows from 
the dispersion relation for Alfven waves u){k) = (Landau and Lifshitz (1987)), and v/ua '> 1 for ultrarelativistic 

particles. 

We shall assume here first the Kolmogorov spectrum normalized at the basic scale fco = 27r/Ao: 

kw{k) = fcou;o(A;/A;o)"^/^ (38) 
where w(k) is the spectral magnetic energy density normalized as fegwo = Bq/Stt. 

We shall adopt from Berezinsky et al. (1990b) (chapter 9) the diffusion coefficient I^n (in the direction of Bq) expressed 
with the frequencies of particle scattering off the waves: 

1/(11, kres) = + V~) = 2lT'^UJBkresi^{kres)/BQ 

D=^-^fdt,^^ (39) 

where v'^ and i/^ correspond to waves propagating along the field and in the opposite direction, respectively. Prom now 
on the subscript || will be omitted. 

Using the formulae above, and performing integration over /j, we obtain 

D{E) = —^vlJ'-^V^\ (40) 



77r(27r)2/3 '^yi 
Numerically 

D(i.) = 2.3xlO"<(15^)"'(^)"'c.V3. (41, 

D{E) from Eq. (41) and from Eq. (35) match together fairly well: at energy Ec they differ by 40%. 

The calculations for other spectra are similar: in case of the Kraichnan spectrum w{k) oc one obtains D{E) oc 

E^/'^ , and for diffusion on shock waves, w{k) (x D = const follows. 

Diffusion coefficient for static magnetic field. 

This case can be considered as scattering off MHD waves in the limit u){k) — > 0. Assuming the spectrum 

with a < 1 for convergence of the integral, we obtain, performing the same calculations as above: 

7r(l - a)[S - a) \ h J 

5.2. Intergalactic medium and formation of turbulent spectrum 

The intergalactic medium is represented by different structures (clusters of galaxies, superclusters, filaments and voids), 
where the medium properties are very much different. In the estimates below we shall normalize all quantities by the 

baryon density n;, = 2.75 x 10~^ cm"'"*, corresponding to 0,1, = 0.044 from WMAP measurements (Spcrgcl et al. (2003)), 
and by the temperature T ~ 10^ K (see e.g. the simulation of Dave et al. (2001) for filaments). Hence, we adopt the 
sound speed 

Cs = {^T/ruHf'^ = 1.2 X lQ\T/l(ff/^cm/s, (43) 

the Alfven velocity 

B , B /2.75 X 10"^ cm-3\^/^ 



and the Coulomb scattering length for electron-electron and proton-proton scattering 

where L ~ 20 is the Coulomb logarithm. The short scattering length Isc in comparison with the basic scale Ic is one of 
the conditions which are needed to provide the turbulent regime. 
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Fig. 3. — Maximum distance to a source rjn&x{E) = 2^ \{E, E'^^^) for diffusive propagation. Continous line: E^'^^ = 1 x 10-^^ eV, dashed 

line: = 1 X 10^2 eV, dotted line: E™^ = 10^3 eV. Magnetic field configuration is {Bo,lc) = (100 nG, 1 Mpc) (left panel) and 

{Bo,lc) = (1 nG, 1 Mpc) (right panel). 



With the characteristics of the media above, we shall discuss now, whether the equilibrium turbulence spectrum in the 
gas can be formed during the age of the universe to- We shall make estimates for the Kolmogorov spectrum. 

The relaxation time to the equilibrium Kolmogorov spectrum is given by Eq. (34) as ~ l/kvk, where Vk is the 
turbulent velocity on the scale k. Taking Vk from pv^ ~ 2kw{k), and using the Kolmogorov spectrum normalized at the 
basic scale as in Eq. (38), we obtain 

1 / \ 1/2 / ; \ -2/3 

1 / p \ I k 



^'^ ko [2koWo) [ko) ■ ^^^^ 



The longest time is needed for the formation of the spectrum at the largest scale fco: 

1 / ^ \ 1/2 



'' = ¥o[2k^o) ■ ^''^ 
Using kowo = Bq/Stt and vq = Bo/\/47rp = ua, we obtain 

TO ~ Ic/UA- (48) 

Numerically it gives 

-0 = 2.4 X 10"^i# ( yr. (49) 



1 Mpc Bo \2.75 X 10-7 

Another estimate of the relaxation time tq can be obtained for the sonic turbulence with the Kolmogorov spectrum. The 
shortest time is given by tq = Ic/Cs^ which for a turbulent plasma is the analogue of the causality condition. Numerically 
it gives 

For a more accurate estimate one can use the relaxation time for the sonic turbulence from Arzimovich and Sagdeev 
(1979) 

Tk ~ . (51) 

kwk ujk 

Using the dispersion relation Uk = k{T/my^'^ and the Kolmogorov spectrum kwk = koWo{k/ko)~'^^^, with koWo ~ nT, 
we obtain numerically an estimate close to that given by Eq. (50). 

From the above estimates we can conclude that the Kolmogorov spectrum cannot be reached for MHD turbulence in the 
voids, where the magnetic field is presumably small, Bq < 1 nG. However, it can be established in filaments (see Eq.(49)) 
in case Bq > 1 nG and /c < 1 Mpc, and it has enough time to be developed in galaxy clusters, where the magnetic field is 
strong and density of gas is larger than in other structures. For other types of turbulence, e.g. for sonic turbulence, these 
conditions can be somewhat relaxed. 



10 




5.3. Some features of diffusive propagation 

Inspired by numerical simulations (Yoshiguchi et al. (2003)) we study the diffusion in various magnetic field configu- 
rations with the basic parameters {Bq, Ic) in the intervals i?o = 10 — 1000 nG and Ic = 1 — W Mpc. As a representative 
configuration wc consider (100 nG, 1 Mpc) with i?c « 1 x 10^° cV. 

The calculated diffuse energy spectra are characterized by the sets {Bq, l^, d), where d is the separation between sources. 
We shall use the following definitions. The case when at all relevant (observed) energies d > ld{E) (or d ^ ld{E)) we 
shall refer to as diffusion in a strong magnetic field. The case when d < ld{E) corresponds to diffusion in a weak magnetic 
field? The extreme case d ^ ld{E) results in the universal spectrum. 
Minimal distance to a source. 

For the diffusion approximation to be valid, a source must be at a distance r > r^i^[E). In principle, velocity of light 
does not enter the diffusion equation, but this equation is not valid when the propagation velocity exceeds the light speed 
c. The value of rmin [E) can be estimated from the condition that the diffusive propagation time must be longer than the 
time of rectilinear propagation, 

2 

r r 
~ D{E) > c' 

as 

r,,,UE) ~ \ld{E). (52) 

For a source at distance r < rram{E) the protons with energy E or higher propagate in a quasi-rectilinear regime. At 
E > Ec one has 

Maximal distance to a source. 

As seen from Eq. (19) the contribution of a source to the fiux at energy E becomes negligible at distances r > rma,x{E) 
with 

rm..iE) = 2^X{E,E^--), (53) 

where E'^^^ is the maximum generation energy provided by a source. 

For the representative case {Bq, l^) = (100 nG, 1 Mpc), rmax is plotted in Fig. 3 as a function of E for different values 
of One can see that in the case of diffusive propagation only nearby sources contribute the UHECR flux. 

In Fig. 4 we have plotted also 

r,s{E) = 2^\{E,Eg/E) 

for different (fixed) ratios Eg/E = 2,5 and 10. The increase of r^s with E provides the increase of diffuse fiux Jp{E) with 
energy when the distance between sources d is large enough (see Fig. 2). 



^ Note, that this definition does not guarantee the (quasi)rectilineaj: propagation of particles. 
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Fig. 5. — Diffusion spectra for two different configurations (Bo, Ic, d). Dashed line corresponds to (100 nG, 1 Mpc, 30 Mpc) with source 
luminosity Lp = 8.5x 10''^ erg/s, while continuous line is shown for the same d but for much stronger magnetic field: (1000 nG, 1 Mpc, 30 Mpc), 
with source luminosity Lp = 1.5 x 10*^ erg/s. In the latter case the GZK cutoff is weak. The spectra are compared with the AGASA data. 



6. DIFFUSE FLUXES IN THE DIFFUSION APPROXIMATION 

In this Section wc will compute the diffuse flux according to Eq. (23) as the sum of fluxes from single sources located in 
the lattice vertexes with a separation d and with a total size a. We will assume the complex generation spectrum given 
by Eq. (8). For the diffusion coefScient we use D{E) = {l/3)cld{E) with ld{E) given by the interpolation formula (29) 
for the case of Bohm diffusion a.t E < Ec- 

At small distances r < r,„in(i?), given by Eq. (52), the fluxes from the individual sources are calculated in the rectilinear 
approximation, and the diffuse flux is given by 

{Eg{E,n))dEg{E,r) 

At large distances r > rinin(-B) the diffuse flux is given by Eq. (23), and with explicit normalization has the form 

where E^ax is the maximum acceleration energy (see Section 4). The calculated spectra are presented in Fig. 5. 

In Fig. 5 the diffuse spectra are shown for the complex generation spectrum (8) and for two sets of {Bo,lc,d) equal 
to (100 nG, 1 Mpc, 30 Mpc), shown by dashed curve, and (1000 nG, 1 Mpc, 30 Mpc), shown by solid curve. Both of 
them are characterized by the same d = 30 Mpc, but in the latter case the magnetic field is much stronger. This case is 
characterized by c? ld{E) at all observable energies, and hence corresponds to diffusion in the strong magnetic field (see 
Section 5). From Fig. 5 one can observe that GZK cutoff is weak in this case. 

The dashed curve in Fig. 5 is characterized by IdiE) = 1{E/Ec)^ Mpc at > E'c = 1 x 10^° eV, and this regime of 
propagation can be described as an intermediate one between that for strong and weak magnetic fields. 

The regime with weak GZK cutoff (propagation in a strong magnetic field) requires a much higher luminosity Lp to 
fit the observational data, Lp = 1.5 x 10*^ erg/s, while the intermediate case (100 nG, 1 Mpc, 30 Mpc) needs only 
Lp = 8.5 X 10^3 gj.g/s_ 

We shall analyze now the regime of propagation in strong magnetic field in more detail. First we will comment 
qualitatively why in the diffusion approximation the GZK cutoff in the spectrmn might be weak or absent. 

The essence of the GZK cutoff consists of much different energy losses above and below 3 x 10^® eV. Consider, for 
example, two protons with energies 1 x 10^^ eV and 1 x 10^° eV, both propagating rectilinearly from a remote single 
source. The first proton looses little energy: the ratio Eg/E (with Eg being the generation energy) is not large, and the 
flux of these protons is weakly suppressed by the generation spectrum. The energy losses of the second proton are large, 
Eg/E is very high and the flux suppression is dramatic (the GZK cutoff). Now consider the diffusive propagation. A 
proton with energy 1 x 10^^ eV, because of the D{E) dependence, travels much longer time than a proton with 1 x 10^" eV, 
and it results in increased ratio Eg/E, making this ratio comparable with that for the second proton. It causes a less 
steep GZK cutoff, or its absence. However, the price for the absence of the GZK cutoff is a very high luminosity of the 
sources Lp, needed to provide the observed flux at > 1 x 10^^ eV, e.g. Lj, > 1 x 10'''^ erg/s for the spectrum shown in 
Fig. 5. 

Let us now come over to the quantitative analysis of absence of the GZK cutoff in the strong magnetic fields. Consider 
the three cases of diffusion regimes at E < Ec with a = 0, 1/3, 1, i.e. D = const, Kolmogorov and Bohm regimes, 
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Fig. 6. — DiflFuse fluxes from lattice-distributed sources with d = 30 Mpc and for magnetic configuration (1000 nG, 1 Mpc). The cases 
a = 0, 1/3, 1 correspond to D = const, Kolmogorov and Bohm diflFusion, respectively at E < Ec ^ 1 x 10^^ eV. Luminosity of a source is 
Lp = 1.5 X 10*^ erg/s. 



respectively (see Eq. (30)). The diffuse spectra calculated for configuration (1000 nG, 1 Mpc, 30 Mpc) are displayed in 
Fig. 6 

Why the three spectra arc so much different at low energies and arc the same at high energies? 

The quantitative explanation can be given explicitly in terms of X{E, Eg), which is the basic parameter of the Syrovatsky 
solution (see Eq. (19)). 

In Fig. 7 we plot the values of 



r^^y,{E,Eg) = 2^X{E,Eg), 

which according to Eq. (19) determines the maximum distance to a source in case the observed energy is E and generation 
energy is Eg. The left panel of Fig. 7 corresponds to D = const, the right panel - to the Bohm diffusion. 

From Fig. 7 one can sec that in the energy interval (1 — 4) X 10^^ eV X{E,Eg) practically does not depend on Eg. 
Then from Eq. (19) one obtains for a single source 



Q{E) 



(7g - l)(47r)3/2 {X{E)f/^l3{E) 1, AX{E) J 
where (3{E) = E-^dE/dt. 

From Eq. (56) and values of A(£') for the two diffusion regimes, one may observe the main effect : strong suppression 
of flux from nearby sources (r ~ d ~ 30 Mpc) in case of the Bohm diffusion and weak suppression in case of D = const 
diffusion. It occurs because the Bohm diffusion coefficient, Db{E) = Do{E / Ec), is small at <C Ec, and hence X{E) 
is also small, which results in exponential suppression of the flux, according to Eq. (56). In case D{E) = Dq, diffusion 
coefficient is large, X{E) is large too (see left panel of Fig. 7), and exponential suppression in Eq. (56) is much smaller. 
This quantitative explanation agrees with the qualitative one, given above: large propagation time, i.e. small diffusion 
coefficient or small A, suppresses the the fiux dX E < Eqzk to the level of the flux at the highest energies. 

Let us come over to the higher energies. 

One can see from Fig. 7 that at > 1 x 10^° eV for D = const diffusion and E >\x 10^^ eV for the Bohm diffusion 
X{E) is small but increases fast with energy up to rinax{E,Eg) ~ 30 Mpc a,t E = i^max- The exponent in Eq. (19) grows 
very fast and all other quantities can be taken out of integral at energy Eg = E'max- After integration one arrives at 
analytic expression 



(56) 



.{r,E) = 



QiErr 



Do ''''^ 4Ao^/^V E 



2_ / -E^max /3(^^max) \ / 



P{E) 



Err 



(57) 



where P{E) = E^^dE/dt and Ao = A(i?, -Emax), which in fact docs not depend on E. From Eq. (57) one can see that 
np{r,E) has universal dependence for all diffusion regimes, provided by D{Ejnax) = Dq{E^^^/ E^Y for all of them. In 
asymptotically high energy regime (not seen in Fig. 6) E^np{E) oc E"^. 

We found and analysed above the diffusion regime with weak GZK cutoff for very strong magnetic field Bo = 1000 nG. 
Can this effect exist for much weaker magnetic field? We have not found such regimes in any realistic cases we studied. 

The strong restriction to existence of these regimes is imposed by upper limit to the distance between sources d. This 
distance cannot be taken arbitrary large. They are limited by maximum acceleration energy i^max, which we keep here 
reasonably high Smax = 1 x 10^^ eV. For rectilinear propagation from a source at distance r the proton with observed 
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Fig. 7. — The values of rmiix(E , Eg) = 2^ \{E , Eg) as function of generation energy Eg for the D = const diffusion (left panel) and the 
Bohm diffusion (right panel). The magnetic field configuration is (1000 nG, 1 Mpc, 30 Mpc). The observed energies E are 1 X 10^® eV - 
curve 1, 3 X lO" eV - curve 2, 5 x IQi'' eV - curve 3, 1 x 10^0 eV - curve 4, 3 x 10^0 eV -curve 5, 1 x lO^i eV - curve 6. 



energy = 1 x 10^" eV imist have at r = 100 Mpc the generation energy Eg = I x 10^^ eV, and thus the sharp cutofT at 
E = 1 X 10^" eV is predicted. For r = 50 Mpc the cutoff energy energy is = 3 x 10^° eV. Assuming r ^ d, one obtains 
d < 50 - 100 Mpc for E^^^ = 1 x 10^^ eV. 

We have performed many calculations with magnetic fields in the range 100 - 300 nG, with Ic in a range 1-10 Mpc, for 
all three difusion regimes a = 1, 1/3, and with in a range 30 - 50 Mpc. In all cases the spectra expose GZK cutoff. 

On this basis we conclude that in case -Emax < 1 x 10^^ eV, i.e. for d < 50 Mpc, the diffusion regime with weak GZK cutoff 
appears only at very strong magnetic fields Bq ~ 1000 nG, and it needs very high source luminosities Lp ~ 10^^ erg/s. 

7. CONCLUSIONS 

We have performed a formal study of the propagation of UHE particles, using an analytic approach. We demonstrated 
that the distance between sources is a crucial parameter which strongly affects the diffuse energy spectrum. 

We have proved that, for a uniform distribution of sources, when the separation between them is much smaller than 
all characteristic propagation lengths, most notably the diffusion length ld{E) and the energy attenuation length latt{E), 
the diffuse spectrum of UHECR has a universal form, independent of the mode of propagation. This statement has the 
status of a theorem and is valid for propagation in strong magnetic fields. The proof is given using particle number 
conservation during propagation and also using the kinetic equation for the propagation of UHE particles. In particular 
the exact solution to the kinetic equation (9) for a homogeneous (i.e. continuous and distance-independent) distribution 
of the sources gives exactly the same spectrum as in the method using particle-number conservation. Note that in Eq.(9) 
the diffusion term is absent due to homogeneous distribution of the sources. 

Another proof of the theorem is given using the diffusion equation (18), and its exact solution (Syrovatskii (1959)) for 
a single source in the case of time-independent energy losses. We calculated the diffuse flux putting the sources at the 
vertexes of a big lattice with size a and with separation d between vertexes (sources). The results of the calculations are 
shown in Fig. 2. One can see that when the distance between sources diminishes from 50 Mpc to 10 Mpc, the spectra 
converge to the universal spectrum (full curves) and at d = 3 Mpc they become identical to the universal spectrum. This 
result is confirmed by analytic calctilations. When the separation between sources is small, stimmation over the sources 
can be replaced by an integration and the corresponding spectrum (21) coincides exactly with the universal spectrum. 

In this paper we have studied the diffusive propagation of UHE particles in intergalactic space, which is considered as a 
turbulent magnetic plasma with baryonic gas density ni, and temperature T, and with two basic turbulence scales Ic, where 
external energy is injected, and ^min, where turbulent energy is dissipated. Turbulent motion is considered to be random 
pulsations at different scales, described by an ensemble of magneto-hydrodynamic (MHD) waves B\exp[i{—u!t + kr)] with 
different frequencies u and different wave numbers k. Diffusion arises due to resonant scattering of particles in magnetic 
fields of MHD waves. At "low" energies, when the Larmor radius is much smaller than the basic scale Ic, the diffusion 
coefficient can be calculated, provided the spectral energy density of the turbulent plasma, w{k), is known as a function 
of k. Such calculations for the Kolmogorov spectrum of turbulence result in the diffusion coefficient given by Eq.(41). 

It is interesting to note, that the diffusion in static magnetic fields can be calculated by this method as the limiting 
case, when wave frequency lu ^ 0. The diffusion coefficient for static magnetic fields is given by Eq.(42). 

The diffusion coefficient in the high energy limit, when rz, 3> Ic, can be reliably calculated as the process of multiple 
scattering. The diffusion coefficient for this extreme case is given by Eq.(35). The diffusion coefficients in the low energy 
regime (r^ "C Ic), given by Eq.(41), and in the high energy regime (rz, ^ Ic), given by Eq.(35), match each other well. 

In practical calculations of the diffuse fluxes we characterize the magnetic configuration by three parameters {Bq, Ic, d), 
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Fig. 8. — Comparison of the analytic diffusive spectrum with the Monte Carlo simulation by Yoshiguchi et al. (2003). The spectra axe 
given in the case of Kolmogorov diffusion with (Bo,'c) = (100 nG, 1 Mpc) and for a single source placed at distance 10 Mpc (left panel) and 
30 Mpc (right panel). 



where Bq is the mean magnetic field on the basic scale Ic and d is the separation of sources. We put the sources at the 
vertexes of a lattice of total size a, for which we typically used a — 300 Mpc. Increasing a docs not change the fluxes. 
The fluxes from the individual sources were found as the Syrovatsky solution (19), and the diffuse flux was found by 
summation over all sources, given by Eq.(23). As a representative case for a strong magnetic field, we considered the 
configuration (Bq) ^c)=(100 nG, 1 Mpc). 

An important feature of the diffusion model is that the observed diffuse flux is produced by nearby sources (see Fig. 
3 and 4). The maximum distance Tmax depends strongly on the maximum acceleration energy i?max- For example, for 
Emax = 1 X 10^^ cV and [Bq, Ic) = (100 nG, 1 Mpc), the maximum distance is less than 70 Mpc, while for i?max = 1 x 10^^ cV 
this distance is 200 Mpc. For smaller magnetic fields these distances are larger (see Fig. 3 right panel). The small radius 
''max of the region, which provides the dominant contribution to the observed diffuse flux of UHEGR, imposes a constraint 
on the diffusion models, which will become more severe when/if particles with higher energies will be observed. As an 
example, let us consider the representative configuration (100 nG, 1 Mpc) and highest observed energy = 3 x 10^° eV. 
In this case we have Ec = ly. 10^° eV, ld{E) 10 Mpc and Zatt = 21 Mpc. To avoid rectilinear propagation from nearby 
sources we should impose the separation d > ld{E) « 10 Mpc. For a maximum generation energy of i?max = 1 x 10^"'^ eV, 
the maximum distance is rmax < 70 Mpc and is only marginally consistent with d > 10 Mpc. The UHECR sources at 
r < 70 Mpc with d > 10 Mpc; and with maximum acceleration energy -Emax ^ 1 x 10^^ cV tentatively imply AGN, whose 
luminosities satisfy the energy requirement for the observed UHECR fluxes in the case of weak magnetic fields. However, 
for the case of a weak GZK cutoff (diffusion in strong magnetic flelds and large separation d) the required luminosities 
are very high Lp ^ 10^^ erg/s. 

The calculated diffuse energy spectra are shown in Fig. 5 for the Bohm diffusion coefficient and for the complex 
generation energy spectrum (8). In case of conflguration with a very strong magnetic gield (1000 nG, 1 Mpc, 30 Mpc) the 
spectrum has a weak GZK cutoff (the solid curve in Fig. 5). This spectrum agrees reasonably with the AGASA excess, 
but requires very large source luminosity Lp = 1.5 x 10^'' erg/s. When magnetic field diminishes to 100 nG the GZK 
cutoff appears (dotted line) . as it should according to the propagation theorem. 

A weak GZK cutoff in the case of diffusive propagation is explained by the fiux suppression at E < Eqzk cutoff due to 
longer propagation time from the source (see Section 6 for the detailed explanation). 

The spectra calculated in this paper in the diffusion approximation are compatible, for the relevant parameters, with 
the numerical simulations in Yoshiguchi et al. (2003). As an example we present in Fig. 8 our spectra from a single source 
calculated in the diffusion approximation for magnetic configuration {Bq, lc)={100 nG, 1 Mpc) with the Kolmogorov 
spectrum of turbulence for a distance to the source of 10 Mpc (left panel) and 30 Mpc (right panel). For the distance 
10 Mpc we also plot the spectrum for the rectilinear propagation shown by the dotted curve. The two curves intersect 
at E'* ~ 3 X 10^*^ eV, as they should provided by the condition ld{E) = r, where r is the distance to the source. At 
E > E^, the particles propagate rectilinearly and the flux is given by the dashed line. This spectrum is compared with the 
muncrical simulations of Yoshiguchi et al. (2003) for the same magnetic field configuration (Yoshiguchi (2004)), shown 
by the crosses. Since the calculations of Yoshiguchi et al. (2003) are not normalized, we equate the fluxes at energies 
with rectilinear propagation as shown in Fig. 8. One can observe considerable disagreement at £^ ~ 8 x 10^^ eV, where it 
is possible to suspect a transitional regime between quasi-rectilinear and diffusive propagation in numerical simulations. 
For the distance to the source r=30 Mpc (right panel) the agreement is much better, probably because at these distances 
the diffusion regime is reached in the numerical simulations. 

The study of this paper is not intended to be a realistic one. We consider the diffusive propagation of UHE particles 
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using an analytic approach with the aim of understanding the basic properties of propagation in strong magnetic fields. 
A realistic propagation should be studied within the hierarchical model of different magnetic fields in different large scale 
structures, as it is done in simulations (Sigl et al. (2003), Sigl et al. (2004), Dolag et al. (2003)). However, we are able 
to reach some practical conclusions. 

Our study is focused on the spectra of UHECR. We demonstrated that the crucial parameter is the source separation 
d. For a wide class of magnetic field configurations, when the separation d is smaller than the diffusion length Id, the 
spectrum is universal, i.e. the same as in case of rectilinear propagation. For the simulation of Dolag et al. (2003) 
the spectrum must be universal. The simulations of Sigl et al. (2003), Sigl et al. (2004) and Yoshiguchi et al. (2003) 
include diffusive and intermediate regimes. According to our calculation, in the diffusion regime very strong magnetic 
fields, B() ~ fOOO nG, and large sc;paration between sources are needed to produce a spectrum with weak (or absent) 
GZK cutoff. In this case very high source luminosities are required, Lp > 1 x 10**^ erg/s. It disfavors diffusion in strong 
magnetic fields as the explanation for the AGASA excess at high energies. 
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APPENDIX 

SPECTRAL ENERGY DENSITY AND SCALE DEPENDENCE OF MAGNETIC FIELDS. 

We shall derive here {B^) oc from the spectral energy density of a turbulent plasma w{k) oc A;"™. The turbulence is 
assumed to be described as an ensemble of hydromagnetic waves with wave numbers k and with vanishing mean electric 

field. The mean magnetic field on each scale A = 27r/fc is {B)\ = B\. i.e. on the scale A the magnetic field is locally 
homogeneous. For the Alven waves the kinetic energy of turbulent fluid is equal to magnetic energy (Landau and Lifshitz 
(1987)). 

Let us write down the Fourier expansion for the wave magnetic field. In the limit ma/c <C 1, where ua is the Alven 
velocity, uj/kuA 1 and the magnetic fields can be considered as quasi-static. Here and below we shall omit the coefficients 
(27r)^, which are inessential for our discussion. 

B{r,t) = j d^kB{k,t)e'^^ 
B*{r,t)= I d^k'B*{k',t)e-'^'^ (Al) 



The energy of magnetic field in the normalizing volume V is 

= ^ / d^'rBir, t)B*{r, t) (A2) 
Putting Eqs.(Al) into Eq.(A2) and using the definition of 5 function 

j d^re'^''-^'^^ = S^{k - k'), (A3) 

we obtain 

W = ^ j d^k\B(k,t)\\ (A4) 



Assuming B{k, t) = B{k, t) we find the energy density of waves and fluid as 

1 

w 

The spectral energy density is then 



w = ^ I dk k'^B'^ik). (A5) 



wik) = ^k'B^k). (A6) 

and hence for the Fourier component B^{k) oc fc^^^, sp = m + 2. Note that the Fourier component B{k) does not have 
meaning of the magnetic field strength over the scale k; moreover the dimension of B{k) is different from a magnetic field. 
The connection between Bx, the average field on the scale A, and the Fourier component B{k), can be readily found from 
relation 



w 



16 



where the integral is to be evaluated over region ^ A. It results in 

4.77- 

Bl ~ ^k^B^{k) = STrkw{k). (A7) 

Using the definition oc /c"**^, we obtain for the physical field B\, sx = m ~ 1, and sp — s\ = 3. 

In particular, for the Kolmogorov spectrum we have sp = 11/3 and sx = 2/3, and thus for defiection of UHE particles 
one should use the wave-number spectrum of magnetic fields in the form 

Bl oc . (A8) 
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